Classical transients and the support of open quantum maps 
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The basic ingredients in a semiclassical theory are the classical invariant objects serving as a 
support for the quantization. Recent studies, mainly obtained on quantum maps, have led to the 
commonly accepted belief that it is the classical repeller - the set of non escaping orbits in the future 
and past evolution - the object that suitably plays this role in open scattering systems. In this paper 
we present numerical evidence warning that this may not always be the case. For this purpose we 
study recently introduced families of tribaker maps [L. Ermann, G.G. Carlo, J.M. Pedrosa, and 
M. Saraceno, Phys. Rev. E 85, 066204 (2012)], which share the same asymptotic properties but 
differ in their short time behavior. We have found that although the eigenvalue distribution of the 
evolution operator of these maps follows the fractal Weyl law prediction, the theory of short periodic 
orbits for open maps fails to describe the resonance eigenfunctions of some of them. This is a strong 
indication that new elements must be included in the semiclassical description of open quantum 
systems. 

PACS numbers: 05.45.Mt, 03.65.Sq 



I. INTRODUCTION 

The strong recent interest in the study of open quan- 
tum systems has evidenced that this area is still in a 
very early stage of development, specially when compared 
with the knowledge that we have of the closed counter- 
parts. As a consequence, it has become an extremely 
active topic in fundamental physics [l|. The most im- 
portant advance concerns scattering systems, where the 
fractal Weyl law (FWL) for the number of long-lived reso- 
nances has been conjectured and tested in different exam- 
ples 0. The corresponding classical invariant distribu- 
tion, i.e. the fractal hyperbolic set of all the trajectories 
non-escaping in the past and future (the repeller), plays 
a fundamental role for the interpretation of the quantum 
spectrum of quasibound (resonant) states. According to 
this law II the number of long-lived resonances scales as 
A rf / 2 , where d is the fractal dimension of the repeller and 
N that of the Hilbert space. After the proof of the frac- 
tal Weyl upper bound for a Hamiltonian flow showing a 
fractal trapped set results that support the validity 
of the conjectured law both for smooth and hard wall po- 
tentials [fj have followed. However, open quantum maps 
provide the ideal arena for this kind of studies, due to 
its computational simplicity @, and accordingly we will 
focus on them in this work. 

Another recent advance in the understanding of scat- 
tering systems has been the development of the short 
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periodic orbits (POs) theory for open quantum maps 
PI- This theory makes use of the shortest POs living in 
the repeller to construct a classically motivated basis in 
which the quantum non unitary operators corresponding 
to the classical maps can be adequately expressed. It has 
been shown that there is a connection between the FWL 
and the possibility to obtain the quantum long-lived res- 
onances by using the detailed classical information of the 
repeller, i.e. the embedded trajectories Q. 

As a result of all this evidence, the repeller is com- 
monly accepted to be the fundamental invariant classical 
structure in phase space able to explain the quantum 
mechanics of open systems. However, in a recent publi- 
cation [9|, it has been shown for a set of maps sharing 
the same repeller and asymptotic features (including the 
escape rate) but having different classical transient be- 
havior, that the quantum properties of the resonances 
can differ strongly from this view. This suggests that 
also transients may play a fundamental role in the de- 
scription of open quantum systems. 

In this paper we further investigate this issue, by pre- 
senting a twofold numerical calculation on the tribaker 
map families introduced in Q. These families differ in 
the way in which the openings arc defined. In the so 
called shift family the area of the opening is the same 
for all of its members, while in the intersection family 
this area changes. First, we have verified that the reso- 
nance distributions follow the FWL for these two kinds 
of maps. Second, we have tested the ability of the theory 
of short POs for open quantum maps to reproduce the 
eigenvalues and the corresponding resonance eigenfunc- 
tions. This approach is used as a tool to investigate the 



actual classical support of the quantum resonances, find- 
ing that the quantum information of some of the shift 
family members cannot be reproduced. These results 
cast doubts on the commonly accepted belief that the re- 
peller is just enough to support the quantization of open 
maps. 

The organization of this paper is as follows: In Sec. 
HI! we introduce the features of the shift and intersection 
families of open tribaker maps. In Sec. IIIII we study the 
effects of classical transients on the spectral statistics. 
The same is done in Sec. IIVI regarding the classical sup- 
port of open quantum maps by using the theory of short 
POs as the main tool. Finally, in Sec. |V] we present a 
final discussion, and our perspectives for future work. 
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II. OPEN MAPS: 



SHIFT AND INTERSECTION 
FAMILIES 



Open maps can be defined as the deterministic evolu- 
tion corresponding to a given (closed) map followed by 
the loss of trajectories which pass through a given region 
in phase-space. The set of trajectories that survive in 
the future and in the past define the backwards and for- 
ward trapped sets, respectively. In turn, the repeller is 
given by their intersection. The structure of this classi- 
cal invariant can be very complex for chaotic dynamics, 
having a fractal dimension in most of the cases. Another 
relevant property of open maps is the decay rate, which 
is defined as the asymptotic loss of probability in one 
iteration. 

In this work we will study the shift and intersection 
families of open tribaker maps defined in that differ 
only in their short time dynamics, sharing the same re- 
peller and decay rate. They are defined preserving parity 
symmetry, and therefore with the same openings in po- 
sition and momentum. 

In ternary notation the map action is given by a 
Bernoulli shift of q — 0. 606162 ■ • ■ and p = 0.e_ie_2£-3 • • • 
(given by the corresponding trits ej = 0, 1,2) as (p\q) = 
. . .e_ 2 e-i.eoeie 2 . . . -t (p'\q') = ■ ■ ■ e_ 2 e_ie .eie 2 . . . 
where the dot is moved one position to the right. Shift 
and intersection families of open maps can be straight- 
forwardly defined replacing the usual trit ej by the open 
trit ij where the value 1 is forbidden (i.e., e, = 0,2). In 
this framework, the shift family members (£>|) are de- 
fined with two open trits corresponding to the fc-th most 
significant trit of both position and momentum. Equiva- 
lently the intersection family members (B l k ) have the first 
k trits open, both in position and momentum. The iter- 
ations of open baker maps in the unit square for k = 1,2 
and both, shift and intersection families are illustrated in 
FiglU 

The quantum version of the shift and intersection map 
families can be defined from the usual quantum tribaker 
map in a Hilbert space spanned by I qutrits (with dimen- 
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FIG. 1. (color online) Iteration of the classical open baker 
map. The unit square divided in 9 regions (3 colors with dif- 
ferent intensities) is shown in the top-left panel. One iteration 
of the colored unit square with the open baker map is shown 
for k = 1 (shift or intersection family) in top-right panel and 
for k = 2 for the shift (bottom-left panel) and intersection 
(bottom-right panel) cases. 
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where the antiperiodic Fourier transform Gjv connects 
position and momentum eigenvectors (\qj) and \pj) with 
j = l,..., AO 
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with N = 1/ (27rfi) (see [HI E| for quantization details) . 

The forbidden 1 in one trit e, can be quantized by 
means of the one qutrit projector ir = I— 11)(1| = |0)(0| + 
|2)(2|. In this way, the projector applied to the i-th qutrit 
can be written as 



with n 
i,j = 1 
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Hi and [IF, IF,] = where 
I. Following this notation, the shift and in- 
tersection families of open quantum tribaker maps can 
be straightforwardly defined as 



b{ = G N n\ 
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with k = 1, . . . , I, and where the form B m i X = HB m i X H 
preserves the parity symmetry of the closed map. 



III. TRANSIENT EFFECTS AND THE 
STATISTICS OF EIGENVALUES OF OPEN 
QUANTUM MAPS 

The geometry of phase space imposes restrictions in 
the asymptotic distributions of the eigcnstates of quan- 
tum systems. In bounded systems of / degrees of free- 
dom, the well known Weyl law establishes that the num- 
ber of eigenstates N(E) of energy less than E is given 
by the number of cells of size h* contained in the corre- 
sponding allowed volume of phase space V(E), i.e. 

N(E) = V(E)/(2Trh) f . (6) 

In open systems the relation between the geometry of 
the phase space and the statistics of eigenvalues is more 
involved. In these systems, the quantum evolution is 
given by a nonunitary operator with right and left de- 
caying nonorthogonal cigcnfunctions. The corresponding 
eigenvalues z n arc complex numbers that fall inside the 
unit circle, that is, \z n \ 2 < 1. The resonances are conjec- 
tured to be localized on the trapped sets fT2l4l4l ] . This 
structure of the phase space has led to the prediction of 
a relation between the number of long-lived resonances 
and the fractal dimension d of the repeller. As previ- 
ously mentioned, the FWL establishes that the number 
of resonances with decay factor fi n = \z n \ = exp[— T n /2] 
scales as ~ N d / 2 . This conjecture has been inves- 
tigated in many systems, including quantum maps, 2-D 
and 3-D billiards, and the modified Henon-Heiles poten- 
tial 
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more scarring has been predicted 
in open systems than in closed ones pjl [ItJ ; this phe- 
nomenon has deep consequences in many applications, 
such as microlasers [l8| . 

The goal of this Section is to test the validity of the 
FWL in the families of open tribakcr maps introduced 
earlier. We remark that all these systems have the same 
repeller, classical information that determines the statis- 
tical behavior of the resonances which is embodied in the 
FWL. In Fig. [5] we show the ordered decay factor /z ra for 
the usual open tribaker map (i.e., the first member of any 

of the two families, B^ 1 ), having Hilbert space dimen- 
sions N — 3 l with I = 4, 5, 6, 7, 8 and 9. In the inset of 
Fig. [2] we show the number of resonances with r n > 0.1 
as a function of N. It is clear that ~ N d ^ 2 with 
d = In (2) / In (3) the fractal dimension of the repeller as 
predicted by the FWL. 

In Fig. [3] we show the same statistics, but for the 
eigenvalues of the open map B%, i.e. the 4th member 
of the shift family. It is worth mentioning that we have 
also studied the other members of this family obtaining 
similar qualitative results as those shown in Fig. [3J We 
clearly see that the ordered decay factor for this map 
using different Hilbert space dimensions roughly behaves 
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FIG. 2. (color online) Ordered decay factor fj, n of the usual 
open tribaker map as a function of n/N. Hilbert dimensions 
N — 3 l with I = 4,5, 6, 7, 8 and 9 (green, blue, turquoise, 
red, gray, and black symbols, respectively). Inset: Fraction 
of states Np with decay rate T > 0.1 as a function of the 
Hilbert space dimension. The solid red line corresponds to 
the prediction of the FWL (N„ ~ N ha(2)/i21n(3)) ). 
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FIG. 3. (color online) Ordered decay factor a n of the open 
map _B| as a function of n/N. The Hilbert space dimension 
is N — 3 l with I = 4,5, 6, 7, 8 and 9 (green, blue, turquoise, 
red, gray, and black symbols, respectively). Inset: Fraction 
of states N^ with decay rate T > 0.1 as a function of Hilbert 
space dimension. The solid red line corresponds to the pre- 
diction of the FWL (JV M ~ ATl«(2)/(21n(3))^ 



as in the case of the usual open tribaker map, i.e. it 
follows the FWL. 

Finally, we consider the statistics of the eigenvalues in 
the case of the intersection family. In Fig. 2] the number 
of resonances with T > 0.1 as a function of N for the 
open maps B % k with k = 1,2,3,4 is shown. We can see 
that again the number of resonances as a function of N 
seems to follow the prediction of the FWL. 

In the following Section we go deeper into the study of 
the classical support of open quantum maps. 
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FIG. 4. (color online) Fraction of states iV M with decay rate 
r > 0.1 as a function of the Hilbert space dimension for the 
open maps B\, with k — 1, 2, 3, 4 (boxes, circles, up and down 
triangles, respectively). 



IV. TRANSIENT EFFECTS AND THE THEORY 
OF SHORT PERIODIC ORBITS FOR OPEN 
QUANTUM MAPS 

The recently developed theory of short POs for open 
quantum maps @, Q (which finds its roots in the cor- 
responding theory for closed systems [l9[) provides the 
ideal tool to test up to which point the repeller, more 
specifically the POs living in it are the actual classical 
support of quantum resonances. The essential ingredient 
of this approach is the open scar function associated to 
each one of these trajectories. 

Taking 7 as a periodic orbit of fundamental period L 
that belongs to an open map, we define coherent states 
\Qj,Pj) associated to each point of the orbit (it has a total 
of L points, all in the repeller) and construct the linear 
combination 

L-l 

= "7f E exp{-27ri(M™ - N8 3 )}\ qj , Pj ). (7) 
VL 3=0 

In this expression m <G {0, . . . , L — 1} and 9j = J2i=o 
where Si is the action acquired by the Ith coherent state 
in one step of the map. The total action is 9l = S 7 and 
A™ = (NS 7 + m)/L. The right and left scar functions for 
the periodic orbit arc defined through the propagation of 
these linear combinations under the open map U (up to 
approximately the system's Ehrcnfcst time r). 

IO = -^jy^^cos (g) (8) 

and 

«J - ^E«i^ e - 2 " A; " tcos (£) • w 

Normalization (N^' L ) is chosen in such a way that 
«J<J = (l#,J<J and (^,J<J - 1. We 



select a number of short POs that approximately cov- 
ers the repeller, whose number scales following the FWL 
prediction Q. Then, we construct an appropriate basis 
in which we can write open evolution operators associ- 
ated to open maps. After solving a generalized eigenvalue 
problem we obtain a classically motivated approximation 
to the system's long-lived resonances Q. 

The short periodic orbit theory is strongly relying on 
the repeller as the source of fundamental information in 
order to describe the long-lived portion of the spectra 
of open maps. In fact, it takes into account the details 
of the open dynamics on the repeller, going beyond the 
consideration of just a measure such as, for instance, the 
fractal dimension which governs the FWL. But now the 
question arises: Is this enough?, or are more elements 
needed? 

We define the performance P of this method 
as the fraction of long-lived eigenvalues that it is 
able to reproduce up to an error given by e = 
v/CMip) - Re(zf °)) 2 + (lm{zf x ) - Im(zf °)) 2 , where 
zf x and zf° are the exact eigenvalues and those given 
by the short POs theory, respectively. Hereafter, we take 
e = 0.001 and consider just the eigenvalues with modu- 
lus greater than 0.01. These threshold values guarantee 
good performance of the approximation and the evalua- 
tion of a meaningful number of eigenvalues in the whole 
range of maps that we have studied. However, if de- 
tailed comparisons are needed, these values should be 
carefully adapted to each map, this going beyond pur- 
poses of the present paper. In Figs. [5] and [6] we show 
P as a function of the number of used POs, Npo s ', here 
our method is applied to the shift and intersection fam- 
ilies for N — 3 5 , respectively. At this relative low value 
of N it is already evident that the usual open tribaker 
map that corresponds to the k = 1 value for both fami- 
lies proves to be well reproduced by this approximation. 
In the case of the shift family the good performance is 
suddenly lost when we arrive at the k = 1/2 threshold, 
i.e., from k = 3 on, the method clearly fails. We notice 
that the resonances of the shift family occupy classically 
forbidden regions for k > 1/2 Q. Nevertheless, the per- 
formance is good for all the members of the intersection 
family. We conjecture that this is because the probability 
outside the region roughly corresponding to the classical 
repeller is erased by means of the projector. To further 
confirm this conjecture we have made the same calcula- 
tion for the I = 6 case, which qualitatively coincides with 
our previous reasoning. The corresponding results can 
be seen in Figs. [7] and [51 

In order to provide further insight on the differences 
among members of these two families of maps and relate 
them to basic features of general open (scattering) quan- 
tum systems we will study the phase space distributions 
associated to the right and left eigenfunctions. For this 
purpose we use a recently introduced representation [2(| 
that reveals the part of the quantum probability which 
is in the closest correspondence to the classical repeller. 
By defining the symmetrical operator hj associated to 
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FIG. 5. (color online) Performance P of the short POs ap- 
proach for the case of the shift family: fraction of long-lived 
resonances found as a function of the number of POs used 
in the calculations Npos- Hilbert space dimension is TV = 3 l 
with I = 5, all possible k are shown. (Color) gray lines with 
squares correspond to k — 1, circles to k — 2, triangles up to 
k — 3, triangles down to k — 4, and diamonds to k — 5. 
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FIG. 7. (color online) Performance P of the short POs ap- 
proach in the case of the shift family: fraction of long-lived 
resonances found as a function of the number of POs used 
in the calculations Npos- Hilbert space dimension is TV = 3 l 
with I = 6, all possible k are shown. (Color) gray lines with 
squares correspond to k = 1, circles to k = 2, triangles up 
to k = 3, triangles down to k = 4, diamonds to k = 5, and 
pentagons to k = 6. 




FIG. 6. (color online) Performance P of the short POs ap- 
proach for the case of the intersection family: fraction of long- 
lived resonances found as a function of the number of POs 
used in the calculations Npos- Hilbert space dimension is 
N = 3 l with I — 5, all possible k are shown. We use the same 
colors and symbols as in Fig. [S] 



the right \Rj) and left (Lj\ eigenstates 

\Rj)(Lj\ 



(Lj\Rj 



(10) 



which is related to the eigenvalue Xj, we construct the 
sum of the first j of these projectors, ordered by decreas- 
ing modulus of the corresponding eigenvalues (\Xj\ ^ 
| Ay | with j <f) 



(11) 



FIG. 8. (color online) Performance P of the short POs ap- 
proach in the case of the intersection family: fraction of long- 
lived resonances found as a function of the number of POs 
used in the calculations Npos- Hilbert space dimension is 
N = 3 l with 1 — 6, all possible k are shown. We use the same 
colors and symbols as in Fig. [7] 



The object to study consists of their phase space rep- 
resentation by means of coherent states |<?,p), which is 
given by 



Qj(q,p) 



\(q,p\hj\q,p}\ 

\(q,p\Qj\q,p)\- 



(12) 
(13) 



i'=i 



In Fig. [9] we show Q32 for the exact resonances and the 
ones given by the short POs approach, obtained for I = 5 
and Npo s = 32 (it is worth mentioning that we have ob- 
served similar results for Qj with j around 32). In the 
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FIG. 9. (color online) Phase space representation of Q32 = 
X]j=i hj f° r the exac t and short POs approach with I = 5 and 
Npos = 32. The map with k = 1 (both shift and intersection 
families) is shown in panels a) for exact eigenstates and b) for 
the short POs approach. For the case of shift family map with 
k — 3 panel c) shows the exact eigenstates and d) the short 
POs approach. The intensities of phase space representation 
are white for zero value and red (dark-gray) for maximum. 



upper panels the case k = 1 (that coincides for both the 
shift and intersection families) makes clear that the short 
POs theory is very accurate in reproducing not only the 
long-lived sector of the spectrum but also the component 
of the resonance eigenfunctions that live on the classi- 
cal repeller. In fact, panels a) and b) are almost indis- 
tinguishable. This shows undoubtedly that the classical 
repeller is the meaningful support of long-lived states in 
this case. Nevertheless, although we do not show here re- 
sults for the intersection family, they all reproduce very 
well these components of the eigenstates, in agreement 
with the conclusions previously shown by means of the 
spectra. When the shift family is analyzed, again a sort 



of phase transition behavior is observed, that spoils the 
short POs theory results in the same fashion as for the 
spectral values. This happens at k > 1/2; here 1/2 = 3, 
case which we show in Figs.|H]c) and d). 



V. CONCLUSIONS 

At present, it is commonly accepted that the classical 
information contained in the repeller is the fundamental 
ingredient to describe quantum chaotic scattering sys- 
tems. However, a recent study has shown that tran- 
sient features of the classical dynamics have crucial effects 
in the resonances of the quantum system. 

In this paper we put this conclusion on firmer grounds. 
First, we have tested the validity of the FWL, one of the 
cornerstones of open chaotic quantum systems theory in 
a family of quantum open tribaker maps @. We have 
clearly shown that the statistics of the eigenvalues fol- 
lows the FWL behavior. In order to unveil the classical 
support of the resonances we have attempted to obtain 
them by means of the short POs approach. This theory 
has been recently linked to the FWL in Ref. Q, and it 
is based on the short POs that live in the repeller. We 
have defined a measure for its performance, showing that 
it fails to describe the quantum resonances associated to 
k > 1/2 for the shift family of maps. In contrast, the 
intersection family, whose opening forces the probability 
to remain roughly inside what can be called the classical 
rcpcller's area, is very well reproduced. 

Summarizing, our results in open quantum maps have 
shown the need to include other ingredients, cither clas- 
sical and/or quantum, additionally to the repeller, in or- 
der to satisfactorily describe the quantum mechanics of 
scattering systems. Further investigations are needed in 
order to find them and extend the short POs theory for 
open quantum maps. 
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